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A stochastic partial differential equation along the lines of the 
Kardar-Parisi-Zhang equation is introduced for the evolution of 
a growing interface in a radial geometry. Regular polygon solutions 
as well as radially symmetric solutions are identified in the de- 
terministic limit. The polygon solutions, of relevance to on-lattice 
Eden growth from a seed in the zero- noise limit, are unstable in 
the continuum in favour of the symmetric solutions. The asymp- 
totic surface width scaling for stochastic radial interface growth is 
investigated through numerical simulations and found to be charac- 
terized by the same scaling exponent as that for stochastic growth 
on a substrate. 



1 Introduction 



The Kardar-Parisi-Zhang (KPZ) equation [1] 

dh d^h A (dhV , , 



is the simplest nonlinear stochastic evolution equation for a growing interface. 
In this equation h[x^ t) denotes the height of an interface at substrate position 
X and time t; ri{x,t) denotes a white noise term and u and A are growth 
parameters related to surface tension and lateral growth respectively. Among 
its applications the KPZ equation models vapor deposition at large length 
scales, Eden growth on a substrate, directed polymers, and one- dimensional 
turbulence. In the case where A = the KPZ equation reduces to the linear 
Edwards-Wilkinson (EW) equation [2] for random deposition. For reviews of 
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the extensive literature on the KPZ, EW and related equations and their 
applicability to evolving interfaces in a wide range of physical phenomena, we 
refer the reader to Refs. [3-6]. 

In this paper we are motivated by one of the above applications, Eden growth, 
to obtain an analogous equation to the KPZ equation in the case when the 
growth is radial from a seed, rather than from the usual substrate. Wc derive 
this equation in Section 3, after recalling the derivation of the KPZ equation 
in Section 2. In Section 4 we peform linear stability analysis for deterministic 
radially symmetric solutions. Numerical solutions arc reported in Section 5, 
along with the exponent f3 governing the width of the interface. We conclude 
with a summary. 



2 KPZ equation for growth on a substrate 



The KPZ equation comprises three terms; surface tension, lateral growth and 
noise. The surface tension acts to smooth out the interface by rounding off 
bumps and filling in hollows. It is assumed that the time rate of change in 
the local height function due to surface tension is proportional to the surface 
curvature 

^ T- (2) 
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The basis of the lateral growth contribution in the KPZ equation is the as- 
sumption that all points on the growing interface move with uniform speed v 
in a direction normal to the interface. The geometry underlying lateral growth 
in the KPZ equation is shown in fig. 1. The lateral growth contribution to the 

time rate of change of the height function is found by projecting the outwards 
normal growth in a direction orthogonal to the substrate. Ignoring overhangs, 
the result which follows from the Pythagorean theorem (fig. 1) is 



(3) 



The noise term in the KPZ equation is assumed to be uncorrelated white 
noise; 

{V{x,t))^0 (4) 

{ri{x, t)ri{x', t')) = 2D5{x - x'){t - t') (5) 
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where D is a, surface diffusion constant. 



The reduction of 
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+ r]{x,t) 



(6) 
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to the KPZ equation, Eq. (1), is acheived by making a small gradient approx- 
imation 1^ <^ 1 and transforming to a co-moving frame. 



3 Continuum equation for radial growth from a seed 

We now consider surface tension, lateral growth and noise in a radial geometry. 
In this case the interface is characterized by its radial position R{9, t) at angle 
9. We assume that R{9,t) is single valued so that there are no 'overhangs' 
along a direction of constant 9. The noise term is now a function of 9 rather 
than X but otherwise it is identical to Eqs. (4), (5). The curvature of the 
surface characterized by R{9, t) is given by 



The geometry for lateral growth in a radial geometry is shown in fig. 2. Simple 
trigonometric manipulations lead to the lateral growth term 




(7) 
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Combining the radial noise term with surface tension proportional to the sur- 
face curvature, Eq. (7), and the lateral growth given by Eq. (8) then gives the 
continuum model for radial growth of an interface: 



The deterministic version of Eq. (9) has radiaUy symmetric solutions R{9, t) — 
f{t) where 

df v , , 



In the absence of surface tension the radially symmetric solutions have the 
form R{t) = vt + C whereas with the inclusion of surface tension R[t) satisfies 
the transcendental equation 

iylog{iy + vR{t)) -vR{t)+vH^C. (11) 



Another class of deterministic 'solutions' are polygons defined by straight hue 
segments 

R{0,t)= f (12) 
cos(6' — c) 



where the parameter c is fixed between the vertices of the polygon. These 
solutions have zero curvature along the faces however they break down at the 
vertices where they are not differentiable. A special class are the regular n 
sided polygons 

^^■^^''^-cos(^-^-^)' —-^-^^ ^^^^ 



with J = 0, .... n — 1. We mention these 'solutions' here because exact regular 
polygons are obtained in on-lattice zero-noise simulations of the Eden A model 
from a seed, where on the square lattice, n = 4 and on the triangular and 
honeycomb lattices n = 6. 



4 Linear stability analysis of deterministic solutions 



The linear stability of solutions R to the continuum model for radial growth 
can be investigated by substituting the perturbed solution R = R + p into 
the growth equation, Eq. (9), and retaining terms to first order in the pertur- 
bation p. This results in a linear partial differential equation to solve for the 
perturbation. In the special case where R — R{t) are the radially symmetric 
solutions the perturbation satisfies 

dp p V d^p 
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Assuming a separable solution p — T(t)H{9) where H{9) — H{9 + 27r) now 
results in 



p{9, t) = ^e"'"W + V (a„ cos(m^) + sin(m^))e(i-"^')'^'"W (15) 



2 m=l 



where 



= / (16) 



is strictly positive. The coefficients a^, &m are determined by the initial con- 
ditions p{9,to). In particular, 



7r 

1^ / P{e,to)^os{m9)d9. (17) 

— TT 

1 ^ 

= ^e(i-m^).a.(to) / ^o) sm(m^)rf^. (18) 

— TT 

Hence if 

TT TT TT 

Jp{9,to)d9= J p{9,to)cos{9)d9 = J p{9,to) sm{9)d9 = (19) 



TTC- 



then the perturbation will decay in time and the radially symmetric solutions 
will be asymptotically stable. 

We can use the above analysis to investigate the linear stability of the regular 
polygon solutions. The idea is to consider an initial regular polygon profile 
as an initially perturbed radially symmetric solution R{t). In this case the 
perturbation initially satisfies 

p{to) ^ Rj{9,to) - R{to) (20) 



where Rj{9,to) is defined in Eq. (13). It is a simple exercise to verify that the 
conditions in Eq. (19) are met in this case and thus the starting regular polygon 
will relax to the radially symmetric solution. This is a direct consequence of 
surface tension suppressing growth at the vertices. The persistence of regular 
polygon solutions in zero noise Eden growth is due to a balance between this 
surface tension and the underlying lattice anisotropy. 
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5 Numerical investigations 

In general the stability of arbitrary profiles can be investigated directly by 
numerically integrating the continuum equation for radial interface growth 
with the given profile as an initial condition. In the numerical studies reported 
below we have employed the simple discretizations 



R(e^ t) = R{iA9, j At) = Rij 


(21) 


7]{e,t)=r]{tAe,jAt)=7]ij 


(22) 


dR Ri,j+1 ~ Ri,j 

dt At 


(23) 


dR R,.+,,, 
89 2A9 


(24) 


d^R _ Ri+ij + Ri-i,j — '^Ri,j 
de^ A9^ 


(25) 



Fig. 3 shows time snapshots (at intervals of 100 time units) of the growing 
interface starting from a diamond in the deterministic limit, r] = 0, with 
growth velocity v — 0.1 for the two cases; {&) u — (sohd line) and (b) 
v — 0.1 (dashed line). In this simulation At = 1/20 and A^ = 27r/500. In each 
case the diamond profile is smoothed out and approaches a circular profile in 
agreement with the predictions of the linear stability analysis above. In the 
integration with surface tension this approach to the circular profile is faster. 

Fig. 4 shows time (number of iterations) snapshots of the growing interface 
starting from a circle in a stochastic integration with rj a random mimbcr in 
the range [—0.008,0.008] and physical parameters v = 1.0, u = 0.1. In this 
simulation At — 1/1000 and A^ = 27r/100. The growing interface in Fig. 4 is 
similar in appearance to the growing interface in Eden growth [7] . 

Fig. 5 shows a log-log plot of the surface width versus time (number of itera- 
tions), averaged over five stochastic integrations, with parameters as in Fig. 4. 
Also shown on this plot is a straight line of best fit with slope P — 1/3. This 
suggests that the continuum equation, Eq. (9), for radial interface growth is 
in the same universality class as the KPZ equation. 



6 Summary 

We have derived the stochastic partial differential equation, Eq. (9), for the 
evolution of a growing interface in a radial geometry. The equation has regular 
polygon solutions as well as radially symmetric solutions in the deterministic 
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limit. Linear stability analysis reveals that the regular polygon solutions are 
unstable in favour of the radially symmetric solutions. The numerical solution 
of the fully stochastic equation indicates that the width of the interface scales 
with the same exponent, (3 — 1/3, as the KPZ equation. Presumably this 
result can be obtained exactly and thus provide explicit confirmation that 
both continuum models of interface growth lie in the same universality class. 
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Fig. 1. Geometry for interface growth from a substrate. The interface is shown as 



the thick grey hne. Here {dh)"^ = (vSt)'^ 
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Fig. 3. Time snapshots of the evolution of an initial diamond shaped interface in the 
deterministic limit of the continuum radial growth equation. The solid line shows 
growth without surface tension and the dashed line is growth with surface tension. 



10 



N = 2'° 



U = 2'^ 




N = 2^*^ N = 2^** N = 2 



Fig. 4. Time snapshots of the evolution of an initial circular interface in the stochas- 
tic continuum radial growth equation. N is the number of iterations with time step 
1/1000. 
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Fig. 5. Log- log plot of the surface width squared versus time (number of iterations) 
in the stochastic continuum radial growth equation. The solid line is the best fit 
line over the domain t G (2.5 x 10^, 10^) with slope /? = 1/3. 
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